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Abstract: Identities involving cyclic sums of terms composed from Jacobi elliptic functions eval- 
uated at p equally shifted points on the real axis were recently found. These identities played a 
crucial role in discovering linear superposition solutions of a large number of important nonlinear 
equations. We derive four master identities, from which the identities discussed earlier are derivable 
as special cases. Master identities are also obtained which lead to cyclic identities with alternating 
signs. We discuss an extension of our results to pure imaginary and complex shifts as well as to the 
ratio of Jacobi theta functions. 
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1 Introduction 



In a recent paper III, we have given many new mathematical identities involving the Jacobi el- 
liptic functions sn(x, m), cn(x, m), dn(x, m), where m is the elliptic modulus parameter (0 < 
m < 1). The functions sn(x, m), cn(x, m), dn(x, m) are doubly periodic functions with periods 
{4K{m),i2K'{m)), {4K{m),i4K'{m)), {2K{m),i4K'{m)), respectively Here, K{m) denotes the 
complete elliptic integral of the first kind, and K'{m) = K{l—m). The m = limit gives K{0) = 7r/2 
and trigonometric functions: sn(a;, 0) = sinx, cn(a;, 0) = cosx, dn(x, 0) = 1. The m —>■ 1 limit gives 
K{1) oo and hyperbolic functions: sn(x, 1) tanhx, cn(a;, 1) —>■ sechx, dn(x, 1) — > sechx. For 
simplicity, from now on we will not explicitly display the modulus parameter m as an argument of 
the Jacobi elliptic functions. 

The cyclic identities discussed in ref. P] play an important role in showing that a kind of linear 
superposition is valid for many nonlinear differential equations of physical interest |31 II]- In all 
identities, the arguments of the Jacobi functions in successive terms are separated by either 2K{m) /p 
or 4K{m)/p, where p is an integer. Each p-point identity of rank r involves a cyclic homogeneous 
polynomial of degree r (in Jacobi elliptic functions with p equally spaced arguments) related to other 
cyclic homogeneous polynomials of degree r — 2 or smaller. In ref. ^j, explicit algebraic proofs were 
given for specific small values of p and r by using standard properties of Jacobi elliptic functions. 
However, identities corresponding to higher values of p and r were only verified numerically using 
advanced mathematical software packages. In this article, we present rigorous mathematical proofs 
valid for arbitrary p and r. As a useful byproduct, we determine explicit forms for the constants 
appearing in various identities. All the identities in ref. corresponded to real shifts of multiples 
of 2K{m)/p or 4K{m)/p. Here, we discuss how to obtain new identities corresponding to pure 
imaginary shifts by multiples of i2K'{m)/p or i4K'{m)/p, as well as identities corresponding to 
complex shifts by multiples of 2[K{m) + iK'{'m)]/p or 4[K{m) + iK'{m)]/p. We also discuss the 
identities for the nine secondary Jacobi elliptic functions like cd (x, m), ns (x, m), ds (x, m). Also, 
we give results for several identities involving Weierstrass elliptic functions and ratios of Jacobi theta 
functions, both of which are intimately related with the Jacobi elliptic functions |2j. In our proofs, 
we classify the identities into four types, each with its own "master identity" which we prove using 
a combination of the Poisson summation formula and the special properties of elliptic functions. 

All our identities involve sums of the following generic form: 

p 

Sp{xo) = J2fi^j) ^ (1) 

i=i 

where /(x) is composed from Jacobi elliptic functions with arguments corresponding to p equally 
spaced points 

Xj = Xo + {j - l)T/p, j = l,...,p , 
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where T is a period of f{x) and the base point xq is an arbitrary complex number. We define the 
quantities P, Q by: 

f{z + 2tK') = {-lff{z), f{z + 2K) = {-lff{z), P,g = 0,1. (2) 

Note that Q = 0,1 correspond to real periods 2K{m),4:K{m) and P = 0,1 correspond to pure 
imaginary periods i2K\m),iAK' {m) respectively. We denote the four possibilities as (+, +), (— , +), 
(+, — ) and (— , — ), where the first sign refers to the sign of (—1)^ and the second to that of (— 1)*^. 
We will derive master identities for each of these four possibilities. 

For example, one of the simplest identities discussed in ref. ^ reads 

p 

^ dn(xj)dn(xj+i) = A , (3) 
i=i 

where A is a constant independent of the base point xq, T = 2K and p is any integer. In this case, 
we have f{z) = dn{z) dn{z + T/p) which corresponds to P = 0, Q = 0, since dn(z + 2K) = dia{z) 
and dn{z + 2iK') = —dn{z). Liouville's theorem can be used to prove the above identity, since dn{z) 
has simple poles within its fundamental region {0,2K,2K + 4iK' ,AiK') at iK' and 3iK' both of 
which we collectively refer to as z* . The identity dn{z* + u) + dn{z* —u) = for arbitrary complex 
u then implies that every pole in the sum is cancelled exactly by a zero of the same order. Thus 
the sum is an analytic function without any poles in the finite part of the complex plane and by 
Liouville's theorem must be a constant This is an explicit illustration of the general principle 
underlying the identities, namely that the orders of poles in a higher order polynomial are reduced 
by some zeros leading to simpler sums. However this method does not yield the constants like A 
explicitly. In fact, using the Poisson summation formula and special properties of Jacobi elliptic 
functions, we show below that the constant A in Eq. (jS)) is given by 

A P r'^A ( \A ( MT/ W (a (OT^I ^ Cn{2K/p)Z{f32K) \ ... 

2K Jo + ^/^) dx = p \ dn{2K/p) sn{2K/p) J ' 

where Z is the Jacobi zeta function {Z = Z{(3g,m)) [2^ with (3g = arcsin(sn(g/p, m)) being the 
Jacobi Amplitude function. 

Identities analogous to Eq. Q also hold for sn and cn. For instance, 

/ N / N P Z{(32k) /p-N 

sn sn Xj+i = 777^^ • (5) 

m sn{2K/p) 

Here T is 2K as this is the periodicity of f{z) = sn{z) sn{z + T/p). The expression as given above 
is valid for all integer values of p > 2, with both sides vanishing when p = 2. 

A further generalization that can be easily treated with the techniques developed below is to 
sums that involve r-th neighbors. We consider the case when r and p are coprime integers and 
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1 < r < p — 1, the other cases being included since identities for any choice of p also include the 
identities for the factors of p. Such a generalization of say Eq. (0} above is 

t d„(.,)dn(.,.,.) = V [M^rKM - . (6) 

Another easy generalization is to identities involving terms consisting of a product of an arbitrary 
number of Jacobi elliptic functions. 

The plan of this article is as follows. In Sec. 2, we derive the master identities which form the 
basis for obtaining all the identities in this paper, and are tabulated in Appendices A and B. Sec. 
3 contains a derivation of identities involving alternating signs. In Sec. 4, we present a collection 
of comments, including some which permit a generalization of all the identities to incorporate pure 
imaginary and complex shifts and to present the identities for Weierstrass functions as well as for 
ratios of Jacobi theta functions. 



2 The Master Identities 

In this section, we derive the four master identities corresponding to Q, P taking on values 0, 1, which 
effectively encompass most of the cyclic identities discussed in ref. [T]. The remaining identities in 
ref. |T] correspond to master identities with alternating signs, and these are discussed in the next 
section. 

For completeness we first derive the finite version of the Poisson summation formula 6] that 
fully exploits the equally spaced nature of the sampling points, and which plays a crucial role in 
subsequent derivations. Since f{x) has a period T, we may expand it in a Fourier series: 

/(^) = ^ Yl a,ke{kx/T) , = f{x)e{-kx/T)dx, (7) 

T k=-oo ^0 

where we have introduced the convenient notation e{x) = exp(27rzx). The required sum may then 
be written as: 

P ]^ oo p 

Spixo) = J2fixj) = Tj; J2 akeikxo/T) ^e(A;j/p). (8) 

j=l k=—oo j=l 



Using the simple identity 



V e{kj/p) = {P ^'^ (9) 
otherwise , 



we get 

Sp{xo) = ake{kxo/T) . (10) 

k, p\k 

Note that we need to evaluate only those Fourier coefficients for which is a multiple of p. 
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2.1 Cases corresponding to Q = 



We first derive the two master identities corresponding to Q = [or equivalently T = 2K], allowing 
P to be either 1 or 0. Consider the rectangle ABCD = {-K, K,K + 2iK', -K + 2iK'). We assume 
that f{z) has a finite number of poles inside ABCD situated at points z^^ = iK' + wT/p, where 
w = 0, ±1, ±2, ... and \w\ < p. Let the principal part of f{z) about the pole z^^ be 



E 



a 



(w) 



1 - ^tuY 



(11) 



which makes this a pole of order L^. 

We now use the fact that f{z) is composed of elliptic functions, and this essentially allows 
evaluation of at- To evaluate a^, for k 0, consider the integral over the rectangle ABCD: 

r~K+2iK' 

f(z)e(—kz/T)dz = ak+ / f(z)e(—kz/T)dz 

JK+2iK' 

= ak + J^^ f{z + 2iK') e{-k[z + 2iK']/T) dz 

au + {-If^'q-^^ f{z) e{-kz/T) dz = a,[l + {-lf+\-^>'] , (12) 



where q = exp{—7iK'/K) is the Jacobian nome The contributions of the vertical segments of 
the integration contour are equal and opposite, and cancel each other. 

On the other hand, the sum of the residues of f{z)e{—kz/T) may also be calculated. The 
residue at the pole is: 

Res [f{z) e{-kz/T)] = Res [f{z) e {-k[z - {iK' + wT/p)]/T) q-^ e{-kw/p) 



q ''Res 




27rtk^Jz-z: 



(13) 



For the second equality, we have made use of the fact that only those for which fc/p is an integer 
need to be evaluated, thanks to the Poisson summation formula Eq. pO|) . 

Define L' = Max{Li, L2, . . . , L^} and 7^ = Yl,^ a^^^ , I = 1, . . . , L' , where we set nonexistent 
a^^^ to be zero. We also set L to be the maximum integer such that 7l is nonzero. If L = 0, there 
are no nonvanishing 7 and the function is regular. Using T = 2K the sum of the residues at all the 
interior poles may be written as: 



Res = q-'' 



1=1 



(/-I) 



2K ) 



(14) 
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Thus 



_ 2TTiq' " ^ 
~ 1 + (-l)^+ig-2fc ^ 

Therefore, Eq. ()10|) now becomes 



ll 



(/-I)! V 2ir 



-2m 



J-i 



(15) 



P 

2K 



7/ 



;et(/-l)! V 2K 



-2m\^-^ 



S r 



2^ 



(16) 



We are in a position to derive two master identities (MI) corresponding to the P = 1 and 
P = cases, which we call MI - I and MI - II respectively. We state for convenience the following 
well-known symmetry properties: 

sn(2 + 27^^) = -sn(2), cn(2 + 27^) = -cn(z), dn(2 + 2/^) = dn(z) , 

sn(2 + 2iK') = sn{z), cn{z + 2iK') = -cn{z), dn{z + 2iK') = -dn{z) . 



2.1.1 MI - I: Case Q = 0, P = 1 

MI - I identities result if there are an odd total number of dn and cn, and an even total number of 
sn and cn functions in f{z), i.e. if one considers terms of the form dn" sn'' cn'^, then a + c is odd 
and 6 + c is even. The primitive function of this type is dn{z) and we consider an "archetypal sum" 
0"! from which identities in this class can be derived: 



X 



(17) 



where p is any (odd or even) integer. We note that in the case of MI - I 

rT=2K 



ao 



f{x) dx = ircyi 



(18) 



as can be seen on integrating f{z) around ABCD, and making use of the antisymmetry about 2iK', 
since we are considering the case P = 1. 

Since dn(2;) has a single simple pole at iK' interior to ABCD with 71 = —i, using the Poisson 
summation formula yields 



k^O,p\k 



f kxo' 

q-^''^ [2K 



(19) 



The above expression for o"i(xo) now allows us to re- write Sp{xo) as given in Eq. (fT^ . yielding our 
first master identity: 



(20) 
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Thus all the sums in this class can be written as sums over the higher order derivatives of the 
function dn(2;). The highest derivative order is one less than the maximum of the orders of the 
function f{z) at all the interior poles. We see that the sums involving the Jacobi functions are 
intimately related to their singularity structure in the complex plane. 
As an illustration, consider the sum 

p 

Sp{xo) = '^dn{xj)dn{xj+i)dn{xj+2) ■ (21) 
i=i 

The relevant function is f{z) = dn{z)dn{z + 2K/p)dn{z + AK/p), with poles at iK', iK' — 2K/p, and 
iK' — AK/p within ABCD. The principal part of the function dn(2;)[dn(2 + 2K / p)dn{z + 4K/p) + 
dn[z — 2K/p)dn{z + 2K/p) + dn{z — 4K/p)dn{z — 2K/p)] around z = iK' determines the 7^. The 
singularity of dn{z) is simple, therefore L = 1. Using the identity dn{z + iK') = —i cs{z) we get that 

71 = —i[cs'^{2K/p) — 2cs{2K /p)cs{4:K / p)] . Substituting this result in Eq. (^0)) gives the identity 

p p 

^dn(xj)dn(xj+i)dn(xj+2) = [cs'^{2K/p)-2cs{2K/p)cs{AK/p)]J2Mxj)- (22) 
j=i j=i 

As another example consider 

p 

Spixo) = ^ sn(xj) cn(a;j) dn(a;j)[dn(xj+i) + dn(xj_i)] . (23) 

The relevant function now is f{z) = sn{z) cn{z) dn{z)[dn{z + 2K/p) + dn(z — 2K/p)\. There are 
three poles, one at iK' and the others at iK' ± 2K/p. To get the quantities 7;, it is convenient to 
consider the principal part of f{z) + f\z + 2K/p) + f\z — 2Klp) around z = iK'. At iK', while 
the product of the three functions sn, cn and dn gives an order three singularity, it is reduced by 
one due to the vanishing of a constant term in the expansion of dn{z + 2K/p) + dn{z — 2K/p) 
around the same point. Thus we get that the maximum order of f{z) is L = 2 while 71 = and 

72 = {—2i/m)ds{2K/p)ns{2K/p) . Substitution in Eq. (j20|l leads to the identity: 

p p 

^ sn(xj) cn(xj) dn(a;j)[dn(xj+i) + dn(a;j_i)] = 2ds(2i^/j9) ns{2K/p) cn(xj) sn(x,) . (24) 
i=i j=i 

Several other identities of this type are given in Appendix A. 
2.1.2 MI - II: Case Q = 0, P = 

This case results when there are an even total number of dn and cn and an even total number of sn 
and cn i.e. if one considers terms of the form dn" sn'' cn^, then both a + c and b + c must be even. 
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In this case the relevant primitive function can be taken as dn^(z) and we define and evaluate the 
following archetypal sum: 

. ^ V- 1 2/ N PE vn^ ^ kq^ ( kxo\ , . 

. g d„ (.,) ^ - + ^ ^ ^- j . (25) 

Here E is the complete elliptic integral of the second kind [21 • We note that in this case 71 is zero, 
as the integral of f{z) around the rectangle ABCD vanishes. 

Substituting Eq. in Eq. (fTHjl . yields the second master identity: 

/ f{x)dx + 272E 

JO 



Thus all MI - II identities have derivatives of (\xi'{z) upto order L — 2. This is also the only 
master identity that has a non- vanishing "constant" (independent of Xq) term on the right hand 
side. The simplest member of this class has already been discussed in Eq. (jH)). We note that the 
relevant function for this identity is f{z) = dn(z)dn(z + 2K/p). There are two poles, one at iK' 
and other at iK' — 2K/p. Thus we can construct dn(2; + iK')[dn{z + iK' + 2K/p) + dn{z + iK' — 
2K/p)] = —cs{z)[cs{z + 2K/p) + c?,{z — 2K/p)] and its principal part around z = Q will give us 
the 7;. The principal part of cs(z) around z = is 1/z. Therefore the only relevant number is 
7i = — [cs(2i^'/p) + cs(— 2i^'/p)] = 0. Anyway we have already observed above that for this class 
7i = from the fact that ABCD is a period parallelogram for f{z). Thus the sum of the principal 
parts cancel and so do all the 7;. Hence using Eq. (j^ . we obtain the identity (H)). In fact we can 
easily generalize using the same argument to a cyclic sum of any even number of dn or sn or cn. 
For instance, 

^dn{xj)dn{xj+r)dn{xj+s)dn{xj+t) = 777? f{x)dx , (27) 

2K Jo 

where f{x) = dn(x)dn(x + r2K/p)dn{x + s2K/p)dia{x + t2K/p). 

As another example, we establish the identity 

p p 

dn2(xj)dn2(xj+i) =AJ2 dn^(a;j) + B . (28) 

Writing f{z) = dn^ (z)dn'^ (z + 2K / p) , there are two poles of order two at iK' and iK' — 2K/p within 
ABCD. We find that L = 2 and 72 = 2cs'^{2K/p). Thus applying the master identity leads to 



(29) 



A = -2cs^{2K/p), B = ^( dn2(t)dn2(t + 2K/p) dt + AEcs^{2K/p)\ . 

2K \jo J 

For an example of this class with L = 3, we prove the identity 

p p 

cn(xj) sn(a;j) [dn^(a:;j+i) + dn^(xj_i)] = A^^ cn(a;j) sn(a;j) dn(xj) . (30) 
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We can derive this using the master identity with f{z) = cn(z) sn(z)[dn^{z+2K/p) + dii'^{z — 2K/p)], 
and we find that L = 3 with 72 = and 73 = (2/m) ds{2K/p) ns{2K/p) . Thus the first derivative 
of dn^(z) will appear in the RHS, which indeed leads to the above identity with the constant 
A = -2ns{2K/p) ds{2K/p). Note that 

r2K i-K 

/ f{t)dt= / [/(t) + /(-t)]cit = 0, (31) 

JQ Jo 

since f(t) is an odd function of t. 



2.2 Cases corresponding to Q = 1 

When Q = 1, the function f{z) has a real period AK. We consider the rectangle ABCD = 
(— e, AK — e, AK — e + 2iK', — e + 2iK'), where e is a small positive number, and integrate around 
this rectangle. Poles occur at iK' + wAK/p and iK' + 2K + wAK/p inside the rectangle ABCD. 
If the principal part around iK' + wAK/p is given by the set of coefficients {7;}, the set around 
iK' + 2K + wAK/p is {-7/}, since f{z + 2K) = -f{z). Also note that 

rAK 

ao = / f{x)dx = , (32) 
Jo 

due to antisymmetry about 2K. Applying the Poisson summation formula and following the same 
procedures as for the previous cases, we get the equivalent of Eq. |T6|l : 

2mp^_jif-2m\'-' A;'-i[l - (-l)fc]g-fc/2 / kxo\ 

We note that S2p{xq) = 0, i.e. the sums in these cases vanish for even values of p. This is however 
a trivial identity since f{xj) = — /(xj+p/2) for j = 1, . . . ,p/2. Thus, for Q = 1, it is sufficient to 
only consider identities where p is odd. 

2.2.1 MI - III: Case Q = 1, P = 

This case applies when f{z) has an even total number of dn and cn and there are an odd total 
number of sn and cn i.e. if one considers terms of the form dn*^ sn^ cn'^, then a + c is even and b + c 
is odd. The relevant primitive function here is sn{z) and the archetypal sum is 

Therefore, using Eq. (jHH) in Eq. (j33|l . we get the third master identity: 

Sp[xq) = Vm}^-^j-—^-^a3{xo) . (35) 
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As an illustration take f{z) = sn'^{z)[sn{z + 4K/p) + sn{z — 4K/p)]. This gives L = 1 with 

7i = (2/m^/^)[ns^(4K/p) — ds{AK/p) cs{4:K/p)] . The resulting identity is 

p p 
J2 sn^ixj) [sn(xj+i) + sn(a;j_i)] = (2/m) [ns^(4if/p) - ds{4:K/p) cs{4:K/p)\ ^ sn{xj) . (36) 
i=i i=i 

An example with L = 2 is provided by f{z) = sn{z) dn{z)[sn{z + 4K/p) cn{z + AK/p) + sn(z — 
4:K/p) cn{z-4:K/p)]. This results in 71 = while 72 = -(2/m^/2)ns(4ir/p)[cs(4fr/p) +ds(4i^'/p)] . 
Therefore the first derivative of sn will appear on the right hand side of the identity which we write 
as: 

p 

sn{xj) dn(xj)[sn(xj+i) cn(xj+i) + sn(xj_i) cn(xj_i)] 

j=i 

p 

= {2/m)ns{AK/p) [cs{AK/p) + ds{AK/p)] J2 cn{xj) dn{xj) . (37) 

i=i 

2.2.2 MI - IV: Case g = 1, P = 1 

This case applies when there are an odd total number of dn and cn and there are an odd total 
number of sn and cn in f{z), i.e. if one considers terms of the form dn"^ sn'' cn'^, then both a + c 
and b + c must be odd. The relevant primitive function here is cn{z) and the archetypal sum is 

Therefore, using Eq. (jHH|l in Eq. (jSHl), we get the fourth and final master identity: 

Sp{xo) = -Tl , ;-i ^4(xo) . (39) 

1=1 - 1)! dxo 

As an illustration, consider f{z) = cn^(2) [cn{z + 4K/p) + cn(2; — 4i^/j9)]. This gives L = 1 with 

7i = (2i/m^^'^)[ds^{AK/p) — ns{AK / p) cs{AK / p)] . The resulting identity is 

p p 

cn^{xj)[cn{xj+i) + cn{xj^i)] = (2/m) ns{4:K/p) cs{4:K/p) - ds^(4/C/p)] ^ cn{xj) . (40) 
i=i i=i 

For f{z) = cn^(z) dn(2;)[sn(z + AK/p) + sn(2; — 4/^/^)], we get L = 2 and 7i = while 72 = 
(— 2z/m^/^) cs(4/r/p) ds(4/r/j9) . The resultant identity therefore involves the first derivative of cn: 

p p 

cn^(xj) dn(xj)[sn(xj+i) + sn(xj_i)] = (2/m) cs(4i^/p)ds(4i^r/p) ^ sn(xj)dn(xj) . (41) 
i=i j=i 

This completes our enumeration of master identities for ordinary sums. Many additional exam- 
ples are given in Appendix A. 
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3 Master Identities With Alternating Signs 



Alternating sums provide an immediate and important extension of the master identities discussed 
above. We consider sums of the form: 

Siixo) = E(-irV(x,) . (42) 

where again f{x) has the properties discussed for the ordinary sums in Eq. We are however 
forced to restrict p to be even in this section and as a resuh we only have MI-I and MI-II master 
identities with alternating signs. One of the consequence of having an alternating sum, as we will 
see below, is the appearance of the simply periodic Jacobi zeta function [7] as an important player. 

To clarify the differences that arise between ordinary and alternating sums we first work out an 
example: 

p 

Sp{xo) = 5^(-l)^-Mn2(a;j)[dn(xj+,) +dn(xj„,)] 

r n ^ 

= 2 [ds{r2K/p) ns{r2K/p) - (-iycs'^{r2K/p)\ E(-l)-'"Mn(xj) . (43) 

Here the spacing is r2K/p. Since r and p are coprimes and p is restricted to be an even integer, 
hence for alternating sums r can only take odd integral values. 

To prove the above identity, in the case r = 1, consider the sums Sp and S~, corresponding to the 
positive and negative signed terms in Sp. We have to take f{z) = dn'^{z)[dn{z + T/p) + dn{z — T/p)] 
with T = 2K. 

Sp{xo) = dn^(a:;i)[dn(a;2) + dn(a;p)] + dn^(a;3)[dn(a;4) +dn(a;2)] H hdn^(a;p_i)[dn(xp) + dn(xp_2)] , 

(44) 

and 

Sp{xo) = dn^(x2)[dn(x3) + dn(xi)] + dn^(x4) [dn(x5) + dn(x3)] H h dn^(xp)[dn(xi) + dn(xp_i)] . 

(45) 

We see that 

Sp{xo) = S^{xq + T/p) 

and 

S^ix^) = E dn'ko + jT/p] [dn(xo + jT/p + T/p) + dn(xo + ]T/p - T/p)] , (46) 

where we have defined p = p/2. The important point to note is that while the above sum appears 
to be in the form of an ordinary sum considered earlier by simply replacing p with p, it is not so, as 
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the function /(x) (which usually depends on p) has remained the same, or equivalently the position 
of the symmetric poles is still at iK' -^T /p, rather than iK' ± T/p. 
Applying the Poisson summation formula we get 



P 

2K 



ao 



k^O,p\k 



kxQ 
2K 



(47) 



Note that we now need for k that is a multiple oi p = p/2 and not merely those that are 
multiples of p. For such k we get upon integrating over the same rectangle ABCD that is relevant 
for type I ordinary identities, 



au = 47i[ns{2K/p)ds{2K/p) - {-l)''/Pcs\2K/p)]j 



2k 



(4J 



This is because at the poles iK' ± 2K/p, f{x) e {kx/2K) has a residue of 2ics^(2i^/p)g~'^(— 1)^/^. 
Now the negative signed sum S~ is related to the positive signed one, by merely a shift in the 
argument by an amount 2K/p. Thus subtracting the two sums leads to a cancellation of the zero 
mode term involving oq and also restricts k/p to be odd integers. We then finally get 



»7r 
2K 



ns{2K/p)ds{2K/p) + cs^{2K/p)] ^ 



k/p=odd 



l + q 



2k 



' kxo' 
2K 



(49) 



A similar evaluation of the archetypal alternating sum can be done with f{z) = dn{z) which is 
simpler as there is only one pole at iK' within ABCD: 



27r ^ 



kxo' 
2K 



(50) 



Therefore the stated alternating identity ()43|) follows. 

We can now generalize these arguments and provide master identities for alternating sums. 
Consider an elliptic function f{z) of real period T satisfying Eq. (1) that has poles at iK' + wT/p 
where w = 0, ±1, ±2, .... and \w\ < p. For both MI-I and MI-II classes we have 



S^ixo) 



- E 

K ^ 



kjp=odd 



kxo' 
2K 



(51) 



thus there are no constant terms, even for type II alternating identities. 

For type I and II identities we can write the a^, the counterpart of Eq. (fT^ as 



2TTiq 



-k L 

^ 1 + (-l)"^+ig-2fc 



li 



{I 



-27Ti 

2K 



i-i 



k 



(52) 
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The difference between the two Eqs. (fT^ and that is crucial, is that 7^ = J2w{~'^)^<^i^^ ■ Thus 
at the pole wT/p, the coefficient of the order / principal part gets weighted by a factor of (—1)'", as 
the residue calculation is restricted to those k where k/p is an odd integer. Therefore for instance 
7i does not in general have the meaning of sum of residues at all the poles. This in turn implies 
that it need not vanish for type II alternating identities. 

Defining the first archetypal alternating sum as in Eq. (j5U|) . we then see that for type I identities: 

Some alternating sum identities of type I are: 

p 

5^(-iy-isn(x,)[cn(x,+i) + cn(x,_i)] = , (54) 



Y^i-iy ^sn{xj) cn{xj) dn{xj)[dn{xj+i) + dn(xj_i)] 
j=i 

p 

= 2ns{2K/p)ds{2K/p) J^i-'^Y'^M^j) cn{xj) . (55) 

i=i 

Turning to the type II alternating identities, it turns out that the primitive function in this case 
is Jacobi zeta function Z{x) rather than dn^(x). To see this, consider the second archetypal sum as 

^2{xo) = ti-^y''Z{x,) . (56) 
On using the Fourier series expansion IHUHI and the Poisson summation formula for Z{x) we get: 

j=i ^ k/p=odd ^ y \ ^-^ / 

On following the steps carried out above then leads to the second master identity: 

Some alternating sum identities of type II are 

p 

^(— l)''^"'^dn(xj)[cn(xj+i) sn(xj+i) + cn(xj_i) sn(xj_i)] 



-{4/m)ds{2K/p) ns{2K/p) ^(-l)^-^Z(a;j) . (59) 
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^(-ly"^ cn%Xj)[cn{xj+i) + cn(xj_i)] = {2/m^)cs{2K/p) ns{2K/p) ^(-l)^"Mn2(a;j) , (60) 
i=i i=i 

p 

^(—1)-'^"'^ cn^(xj) sn(xj)dn(xj)[cn(xj+i) + cn(xj_i)] 
i=i 

= -(4/m2)ds2(2ir/p) cs{2K/p) ns{2K/p) ^(-l)^-^Z(xj) 

i=i 

p 

+ (2/m) cs{2K/p) ns{2K/p) ^(-l)^"^ cn(xj) sn(xj)dn(xj) . (61) 

i=i 

Summarizing, for functions of the form f{z) = h{z)[g{z + T/p) + (7(2; — T/p)], which occur 
in ordinary sums, we may use the symmetrized form h{z)[g{z + T/p) + g{z — T/p)] + g{z)[h{z + 
T/p) + /i(2; — T/p)] and evaluate its principal part at iK'. On the other hand, for alternating sums, 
we may use the anti-symmetrized form h{z)[g{z + T/p) + g{z — T/p)] — g{z)[h{z + T/p) + /i(2; — 
T/p)] and consider its principal part at iK'. Its generalization to more complex forms of f{z) 
is straightforward. Using the master identities derived in this and the previous section and this 
methodology, we have obtained a large number of identities, some of which are given in Appendices 
A and B. 



4 Comments and Discussion 

In this section, we give some general comments and extensions in several new directions. 

(i) Identities for Auxiliary Functions: Until now, we have discussed identities for the three 
basic Jacobi elliptic functions sn, cn, dn. However, nine auxiliary functions are also frequently found 
in the literature. They are: 

, 1 , cnu snu 

nau = ; cau = ; sau = ; 

dn M dn M dn m 

1 CUM dnu 

nsu = ; csu = ; dsu = ; 

sn M sn M sn M 

1 , dnu snu 

ncu = ; dc M = ; sc m = . (62 1 

cnu cnu cnu 

Identities for these auxiliary functions are readily obtained via the following relations j2l |H] : 



dn(M, m) = y/l — mnd {u — KjUi) = —i cs {u — iK', m) = z-\/l — msc {u — K — iK' , m) , (63) 
sn(u, m) = cd (m — K, m) = —= ns (u — iK' , m) = —= dc(u — K — iK' , m) , (64) 
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cn(M,m) = — Vl — msd {u — K,m) = —= ds (u — iK' , m) = = nc (u — K — iK' ,m) . (65) 

wm win 



As an example, consider the identity dn(M,m)dn(u + K,m) = — m. Using Eq. (jUHjl . we 
obtain 

nd (x, m) nd (x + K,m) = —^^=^= , 

V 1 — m 

cs (x, m) cs (x + i^, m) = — yl^^m , 

sc (x, m) sc (x + i^T, m) = . (66) 

V 1 — m 

(ii) Identities for Pure Imaginary Shifts: So far, we have focused our attention on identities 
involving Jacobi elliptic functions evaluated at points separated by real gaps T/p, with real T. As 
mentioned in ref. PP, since Jacobi functions are doubly periodic, we can convert each identity to 
another one involving points separated by pure imaginary gaps iT' /p, with real T'. The procedure 
consists of taking any given identity, writing it for modulus 1 — m [noting that K{1 —m)= K'{m)], 
using the standard results [3 E] 

sn(x, 1 — m) = dia{ix + K{m) + iK' {m),m) , 

^/l — m 



1-y/fYl 

cn(x, 1 — m) = cn{ix + K{m) + iK'{m),m) 

yl — m 



dn(x, 1 — m) = ^/msn{^x + K{rn) + iK'{rn),rn) , (67) 

and changing to a new variable u = ix + K{m) + iK'{m). 

For instance, again consider the simple identity dn(x, m) dn(x + K^rn) = — m. Re-writing 
with modulus 1 — m, using Eq. (jU7j) . and changing to the new variable u = ix + K{m) + iK'{m) 
gives a simple identity involving a pure imaginary shift 



sn(M, m) sn(u + zfr'(m), m) = l/i/m • (68) 

A more non-trivial example consists of identity (42) in ref. P]: 

sn(x, m) sn(x + AK{'m)/3, m) sn(x + 8fr(m)/3, m) 

= [sn(x, m) + sn(x + 4_ft'(m)/3, m) + sn(x + 8ii'(m)/3, m)] , (69) 

where q = dn(2JC(m)/3, m). 
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The corresponding identity with pure imaginary shifts is 

dn{u, m) dn{u + AiK'{m) /3, m) dn(u + 8iK'{m) /3, m) 

— (1 — TTl) 

= ^ [dn(M, m) + dn(M + AiK'{m) /3, m) + dn(M + 8iK'{m) /3, m)] , (70) 

where q' = dn{2K'{m)/3, 1 — m). 

(iii) Identities for Complex Shifts: Just as we have derived identities containing pure imaginary 
shifts, we can also derive new identities involving complex shifts. Here, the procedure consists 
of taking any given identity for real shifts, writing it for modulus 1/m (noting that K{l/m) = 
^/rri[K{m) + iK'{m)] [2]), using the standard results 

sn(x, — ) = y/m sn(— m) ; cn(x, — ) = dn(— m) ; dn(x, — ) = cn(— m) , (71) 
m \/m m \ m m \ m 



and changing to a new variable u = xj ^fm. 

As a simple example, let us once more take the simple identity dn(x, m) dn(x + -ft', m) = \J\ — m. 
It now transforms to 



cn(M, m) cn(u + K[m) + iK'{m), m) = —i ^ 'm^ ^ ^^^-^ 



As a second example, take identity (45) in ref. P: 

cn(a;, m) sn(x + 4:K/3, m) sn(x + 8K/3, m) + cn(a; + 4:K/3, m) sn{x + 8K/3, m) sn(x, m) 
+cn(x + 8K/3, m) sn(x, m) sn(x + 4K/3, m) 

= ^ [cn(a;, m) + cn{x + AK/3, m) + cn{x + 8K/3, m)] , (73) 

m 

where q = dn{2K{m)/3,m). 

The corresponding identity with complex shifts is 

dn(ii, m) sn{u + A{K + iK') /3, m) sn(M + 8{K + iK') /3, m) 
+ dn(M + 4(K + iK') /3, m) sn(M + 8{K + iK') /3, m) sn(M, m) 
+ dn(M + 8(K + iK') /3, m) sn(M, m) sn(M + 4(ir + iiT') /3, m) 

= -(1 + r)2 [dn(M, m) + dn(u + 4(is: + iK')/3, m) + dn(u + 8{K + zis:')/3, rn)] , (74) 
where r = g(l/m) = cn(2{i^'(m) + zi^r'(m)}/3, m). 

(iv) Identities for Ratios of Jacobi Elliptic Functions: In applications involving linear su- 
perposition of solutions of nonhnear differential equations [Sllll, one often needs identities for ratios 
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of Jacobi elliptic functions like cndn/sn. These can be obtained from the identities derived in this 
paper. For example, noting that 

cnx dnx dn2x + cn2x 



snx 



sn 2x 



m cn(2x + iK', m) + dn(2x + iK', m) 



(75) 



gives rise to 



cn {x + 2K/3) dn (x + 2K/3) cn (x + AK/S) dn (x + AK/S) 



sn {x + 2K/3) 



sn {x + 4:K/3) 



cn{x + 4K/3) dn{x + 4K/3) cnx dnx ^ cnx dnx cn (x + 2/^/3) dn (x + 2fC/3) 



sn(x + 4i^'/3) snx ' snx sn{x + 2K/3) 

= -m[cn(2M) cn(2u + AK/3) + cn{2u + AK/3) cn(2u + 8K/3) + cn{2u + 8K/3) cn(2u)] 

-[dn(2M) dn(2M + 4K/3) + dn(2M + AK/3) dn{2u + 8K/3) + dn{2u + 8K/3) dn{2u)] , (76) 

where u = x + iK'{m)/2. In the above derivation, the cndn terms cancel in view of identity (33) in 
ref. [Tl. Further, the right hand side of Eq. (f7H|l has the constant value g(2 + g) [m — (1 + g)^]/(l + 
qY , q = dn{2K{m)/3,m), due to identities (32) in ref. p. 

Other identities involving ratios follow from useful equations analogous to Eq. ()75|) : 



snx dnx 

cnx 
snx 



cn2x snx cnx 1 — dn2x 



cnx 



sn2x dnx msn2x snx dnx 

dn 2x — cn 2x dn x 1 + dn 2x 



1 + cn 2x 
sn 2x 



cnx dnx (1 — m)sn2x snx cnx sn2x 



(77) 



(v) Cyclic Identities for Weierstrass Functions: Jacobi elliptic functions are closely related 
with the Weierstrass function V{u) OEI, the relations being: 



snu = [V{u) — 63 



cnu 



V{u)-ei 
V{u) - 63 



1/2 



dnw 



V{u)-e2 
V{u) - 63 



1/2 



(78) 



where 



ei = (2 - m)/3 , 62 = (2m - l)/3 , 63 = -(1 + m)/3 . (79) 

'P{u) has implicit arguments corresponding to its two periods 2uji = 2K[m) and 2ci;3 = 2iK'[m) 
[7j. Using this relationship and identities obtained by us, we can immediately write down identities 
for Weierstrass function for shifts in the units of cui/p, uj^/p and UJ2/P where UJ2 = —{uJi + cua) . For 
example, using identity (I2H1), one gets 

j2V{u + 2{j-l)u,/p)V{u + 2ju,/p) = {B + pAe,-pel)-{A-2ei)j2'Piu + 2{j-l)ij,/p) (80) 



where A, B are the constants appearing in Eq. 
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(vi) Cyclic Identities for Jacobi Theta Functions: The connection between the four Jacobi 
theta functions 6i{z),92{z), 63(2), 94(2) and the Jacobi elliptic functions is given by [7j 

snu = — — - ^ , , , cnu = ^ , , , an u = (1 — m) ' ^ , , , 81 

where z = 2K{m) ■ Therefore, any of our cyclic identities for real, imaginary or complex shift can 
also be re-written as identities for the ratios of Jacobi elliptic functions for shifts in units of Tc/p, 
7tt/p or 7r(l + t)/p respectively where r = iK' / K. As an illustration, we consider identity ()85|) . In 
terms of theta functions, one gets 

^ e,{z + {j-l)^lp) ^ el{2nK/p) \ ^ es{z + {j-l)7r/p) 

fAe,{z + {j-l)7r/p) I i\ e!{2nK/p) ^,e,{z + {j - l)7r/p) ' ^ > 



It is a pleasure to thank Profs. Henry McKean and Ashok Raina for insightful comments and 
suggestions. The work of AK and US was supported in part by grant FGO2-84ER40173 from the 
U.S. Department of Energy. 
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Appendix A: Examples Using the Master Identities 



In this Appendix, we present a collection of identities involving cyclic combinations of Jacobi 
elliptic functions. These identities are derived by using various choices of f{x) in the four master 
identities [Eqs. (pUjl . (j^ . (jH^j) . (jH^ ] developed in the text. We use the notation a = r2K/p and 
b = rAK/p, where 1 < r < p — 1 and {r,p) = 1. We also use a' = s2K/p, a" = t2K/p and 
h' = s4:K/p, where s,r,t are all distinct. Also we use Sj = sn(xj) etc., where Xj = xq + (j — l)T/p. 
Note that T = 2K for the first two master identities, while it is AK for the remaining two. We note 
that these identities are not exhaustive (as indeed they cannot be) but are meant to be representative 
low L identities. 

Examples belonging to the class MI - I 

L = : 



i=i 



(83) 



L = 1 



\l-l)/2 {l-l)/2 I 

n cs2(A;a) + 2(-l)('-^)/2 ^ ^ cs{[n - k]c 

k=l k=l n=l,n^k 



which is valid for any odd integral / < p. This is a generalization of Eq. ()22|) which corresponds to 
/ = 3. In the special case oil = p this identity takes the simpler form 



j=l n=l P j=l 

A special case of this identity for p = 3 has been worked out in PJ. 



(85) 



E d'^{dj+r + dj^r) = 2[ds(a) ns(a) — cs^(a)] E dj ■ 
j=i j=i 

p p 
E Cj{cj+rdj+r + Cj-rdj-r) = (— 2/m)cs(a) [ds(a) — ns(a)] ^ dj . 



(86) 
(87) 



1^— 2/m)cs(a)[ds(a) — ns(a)] dj 
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p p 
dj{dj+rdj+s + dj^rdj-s) — —2 [cs(a)cs(a') + cs(a — a'){cs{a) — cs(a')} ] ^ 



o?j(cj+rCj+s + Cj-r-Cj-s) = (— 2/m) [ds(a)ds(a') + ds(a — a'){cs(a) — cs(a')}] ^ o?j 



p p 
+ Sj-rSj-s) — (2/?7^) [ns(a) ns(a') + ns(a — a'){cs{a) — cs(a')}] ^ 

j=l 3=1 



Cj(cj+r(^j+s + Cj-rdj-s) — (—2/m) [{cs(a') + cs(a — a')}ds(a) — ds(a — a')ds(a')] ^ (ij 



p p 
y^ Sj(sj+r(ij_i_s + Sj-rdj-s) = (2/m) [{cs(a') + cs(a — a')}ns(a) — ns(a — a')ns(a')] ^ (ij 



'^s'^[cj^rCj+sdj+t + Cj_rCj-sdj-t] = (2/m^) cs(a)cs(a")ds(a') ns(a) 

+cs(a')cs(a")ds(a) ns(a') + ds(a)ds(a')ds(a") ns(a") - cs(a - a")ds{a - a') ns^(a) 

p 



—cs{a" — a') ds(a — a') ns (a') — ds(a — a")ds(a' — a") ns (a') 



p p 

y^ (i^[cj+j.Sj+j. + Cj_rSj_r] = — 2[cs^(a) + ds(a) ns(a)] ^ CjSj 

p p 
y^ CjSjdj[dj+r + c^j-r] = 2ds(a) ns(a) ^ CjSj . 
j=i j=i 

p p 
y^ Sjdj{cj^rdj+r + Cj^rdj-r) — — 2cs(a) [ns(a) + ds(a)] ^ CjSj 
j=i j=i 

p p 
J2 + = -(2/m)cs(a)ns(a) J2 CjSj ■ 

p p 
E«j(c|+^ + 4-^) = (2/m)cs(a)ds(a) ^CjSj . 
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p p 

SjCj{dj+rdj+s + dj^rdj-s) — — 2cs(a)cs(a') ^ CjSj . (100) 

j=l 3=1 

P P 

SjCj{cj+rCj+s + Cj-rCj-s) = (— 2/m)ds(a)ds(a') ^ CjSj . (101) 

j=l 3=1 
P P 

SjCjisj+rSj+s + Sj-rSj-s) = (2/m) ns(a) ns(a') ^ c^-Sj . (102) 

3=1 3=1 

p p 

djCj{sj-^rdj+s + Sj-rdj-s) — —2 ns(a) cs(a') ^ CjSj . (103) 

p p 

djSj{cj^rdj+s + Cj-rdj-s) = "2 ds(a) cs(a') ^ c^Sj . (104) 

o?^(sj+rCj+s + Sj-rCj-s) — 2 [cs(a) ds(a — a') — cs(a')ns(a — a')] ^ CjSj . (105) 

3=1 3=1 



dj'^j^j{dj+r "I" ^j-r) 

cs^(a) ns^(a) + ns^(a)ds^(a) + ds^(a)cs^(a) + 3cs^(a) ns(a)ds(a)J ^ c^Sj . (106) 



3=1 

= -2 



3=1 



L = 3 : 



4('^j+r- + dj_r) = 2 ns(a)ds(a) ^ ^ + 2cs^(a)[cs2(a) - ns(a)ds(a)] ^ ci^- . (107) 

3=1 3=1 3=1 



E4(^.V + 4-.) = -2cs^(«)E4 

i=i 3=1 



+2 [cs2(a) ns2(a) + ns2(a)ds2(a) + ds^{a)cs'^{a) - 3cs^{a) ns(a)ds(a)J ^ dj ■ (108) 



y^ CjSjdj{cj+rSj+r + Cj-rSj-r) — (2/m^)ds(a) ns(a) E "I" (2/"^) 



cs (a) ns (a) 



+ns^(a)ds^(a) + ds^(a)cs\a) - ns(a) ds(a)(cs^(a) + ns\a) + ds\a)) 



J2d,. (109) 

3=1 
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L = 4 



p 



'^d'^{sj+rCj+r + Sj-rCj-r) = —2 ns(a)ds(a) ^ SjCj(i| + 2cs^ (o) [cs^(a) + 3 ns(a)ds(a)] ^ SjCj . (110) 
j=i j=i j=i 



p 

■2 



Yl dj{sj+rCj+s + Sj^rCj-s) = "2 ns(a)ds(a ) ^ SjCjdj 

r 1 ^ 

+2 [cs(a)ds(a)cs(a') ns(a') + ns(a)ds(a'){cs^(a) + cs^(a')}J X] ^j^j ■ (m) 

Examples belonging to the class MI - II 

L = : 

p 

J2djdj+r = P (dn(a) - cs{a)Z{p2rK)) ■ (112) 



i=i 

p 



SjSj+r = — Cs{a)Z{l32rK) ■ (113) 
j=l ^ 

^ / ds(a)Z(^\ 

Cj9+r = P cn(a) . (114) 

3=1 \ m J 

It may be noted that the last two identities are vahd for p > 2 while for p = 2 both sides vanish 
identically. Further, the product of any even number (< p ) of dn's, sn's or en's is also a constant, 
the constant being the integral of the corresponding function over the period 2K. For example. 



P p i'2K 

djdj+i ■ ■ ■ dj+r-1 — / dn(a;)dn(a; + 2K/p) ■ ■ ■ dn{x + (r — l)2K/p) dx, for r even . (115) 
IK Jo 



For dn's of course even the product of all p of them (i.e. did2...dp) is a constant =(1 — m)^/^ in case 
p is even. 



J2 (^jdj{Sj+r + Sj-r) = . (116) 



.7=1 



H djSj {Cj+r + Cj-r) = . (117) 
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p 

CjSj{dj+r + dj-r) = . (118) 

3=1 

P 

Cj{dj+rSj+s + dj_rSj-s) = . (119) 

P 

Yl (dj+rCj+s + dj^rCj-s) = . (120) 

i=i 

p 

5Z dji^j+r^j+s + ^j-rSj—s) = . (121) 

i=i 



L = 2 : 



J2dy,+r = Aj2d] + B, (122) 
where ^ = -2cs2(a), B = ^ (/^^^ dn\t)dn^{t + 2rii'/p) (it + 4Scs2(a)). 

J2cjSj{Cj+rSj+r + C^-r-Sj-r) = I] ^ + ^ I / f {x)dx + 272 E j , (123) 

j=l j=l \ ^ / 

where f{x) — cn{x) sn(x)[cn(x + a) sn{x + a) + cn{x — a) sn{x — a)] and 72 = — (4/m^) ns(a) ds(a). 

5^ Cjdj{cj+rdj+r + Cj^rdj-r) ^ c?? + ^ (^y^ /(2;)(ia; + 272£;j , (124) 

where f{x) — cn(a;)dn(a;)[cn(a; + a)dn{x + a) + cn{x — a)dn{x — a)] and 72 = (4/m)cs(a)ds(a) . 

where f{x) = sn(a;)dn(a;)[sn(a; + a)dn(a; + a) + sn(a; — a)dn(a; — a)] and 72 = (— 4/m)cs(a) ns(a) . 

E d]{d]^rdj+2r + d]_,d^_2r) = '12 E ^5 + ^ X /(^)^^ + ^72^ j , (126) 

where fix) = dn^(a;)[dn^(a; + a)dn(,r + 2a) + dn^(x' — a)dn(x — 2a)] and 72 = 2cs(a)[cs'^(a) + 
2cs(2a)ds(a) ns(a) + cs(a)ds(2a) ns(2a)] . 

P P „ / r2K \ 

Yd]{d,^r + d,.r) = -l2Y.d] + i^yj^ f{x)dx+2-f,Ej , (127) 
where f{x) — dn^(x)[dn(x + a) + dn(x — a)] and 72 = — 2ns(a)ds(a) . 
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i=i i=i ^-^ V-^o 

where f{x) = sn^ {x)[sn{x + a) + sn{x — a)] and 72 = — 2cs(a)ds(a)/m^ . 

E + 9-,) = -72 E ^? + ^ / /(^)^^ + 272^ 

j=i i=i ^-f^ 

where f{x) = cia^{x)[cn{x + a) + cn{x — a)] and 72 = — 2ns(a)cs(a)/m^ . 

E + = -72 E ^? + ^ / /(^)^^ + 272i^ 

j=l j=i 
where /(x) = dn^(x)[dn^(a; + a) + dn^{x — a)] and 72 = 12cs^(a) ns(a)ds(a) . 

E + = -72 E ^? + ^ f + 272^ 

i=i i=i 

where f{x) — sn^(x)[sn^(x + a) + sn^(x — a)] and 72 = — 12ns^(a)ds(a)cs(a)/m' 

E + cj-.) = -72 E ^? + / /(^)^^ + 272^^ 

i=i i=i ^-^ V-^o 

where f{x) = cn'^(a;)[cn^(a; + a) + cn^(a; — a)] and 72 = 12ds^(a)cs(a) ns(a)/m^ . 



(128) 



(129) 



(130) 



(131) 



(132) 



P P p / i'2K \ 

E CjSjdj{Cj+rSj+rdj+r + Cj^rSj-rdj-r) ^ E + ^ (^/q f{x)dx + 2-f2E j , (133) 

where f{x) = cn(a;)sn(a;)dn(x)[cn(x + a)sn(x + a)dn(x + a) + cn{x — a)sn{x — a)dn(x — a)] and 
72 = {-4/m^)[ns\a){cs\a) + ds\a)) + cs\a)ds\a)] . 

L = 3 : 



E (^jSjdj{d^j+r + ^ -2cs2(a) E CjSjdj 



i=l 



E (^j^'jdj{sj+r + Sj-r) = -(2/m) cs(a)ds(a) E cjSjdj 



p p 
E Cj(ijSj((ij+r + c^i-r) = — 2ns(a)ds(a) '^CjSjdj . 

3=1 3=1 

P V 

E SjC^jdjicj+r + Cj-r) = (2/m) cs(a) ns(a) E CjSjC^i 



(134) 
(135) 
(136) 
(137) 
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E «iCj4(4+r + = -4 cs'(a) ns(a)ds(a) £ c,s,d,- 



-4/m^) ns^(a)ds(a)cs(a) '^CjSjdj . 



p 



p 



E ^j-«jS^(c|+r + = (-4/m2) ds\a)cs{a) ns(a) ^ CjSjfij 



(138) 

(139) 
(140) 



r 1 
E CjSj<ij((ij+r + c^j-r) = 2 cs'^(a) — ns^(a)ds^(a) — ds^(a)cs^(a) — cs^(a)ns^(a) E'^i'^J^i • i^^^) 



L = 4 : 



i=i 



2m 

m-t2 - — (m + 1)74 



2K 



p2K 

/ /(x)cix + 272^; 
Jo 



+ p I -72 + 



m74\ 
3 J ' 



(142) 



where 72 = — (l/3m^) cs(a)ds(a) 5 + 5m + cs^(a) + ds^(a) + 4ns^(a) , 74 = — (2/m^) cs(a)ds(a) 
and /(x) = sn^{x)[sn{x + a) + sn(a; — a)]. 



j-r) 



2 \ ^ 4 I 

^ 74 E + 

r2K 



2777- 

"^72 — ^("^ + 1)74 



E' 



p 



2K 







f{x)dx + 2^2E 



+ P(-72 + ^) , (143) 



where f{x) = dn^{x)[dv?{x + a) + dii^(a: - a)], 74 = -2cs^(a) and 72 = -(4/3) cs^(a)(m - 2) 
2[cs^(a) + cs2(a)ds^(a) + ds^(a) ns^(a) + ns2(a)cs2(a)]. 

Examples belonging to the class MI - III 

L = : 



(144) 



L = l : 
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There are identities in sn and cn analogous to the one given by Eq. ()22|) and its generahzations 
in Eq. ()84|) and Eq. ()85p . In particular, for odd I < p we have: 



p 



J2sjSj+r---Sj+(l-i)r = (145) 

(/-l)/2 {l-l)/2 I \ p 

_l)a-i)/2 -Q ns2(A;6) + 2 J2 n ns([n-A;]6) j^'^r 

k=l k=l n=l,n^k / j=l 



m 



{l-l)/2 



When / = p, the resulting identity takes the simpler form 



^(p-l)/2 p 



i=l \ n=l P J j=l 

X: s'^{s,+r + s,.r) = (2/m)[ns2(6) - ds(6)cs(6)] s, . (147) 
i=i i=i 

^ Cj(cj+rSj+r + Cj_rSj_r) = (2/m) ns(6)[cs(6) — ds(6)] ^ Sj . (148) 

p p 

dj{dj+rSj+r + dj^rSj~r) = 2 ns(6) [— cs(6) + ds(6)] ^ Sj . (149) 
i=i i=i 

p p 
^ Sj(cj+^,Cj+, + Cj_rCj-s) = (-2/m) [ds(6)ds(6') + ds(6 - b'){ns{b) - ns(6')}] E • (^^O) 
j=i j=i 



p p 

Sj{dj-^rdj+s + dj^rdj-s) = —2 [cs(6)cs(6') + cs(6 — 6'){ns(6) — ns(6')}] ^ Sj . (151) 



p p 

Sj(sj_|_r'Sj+s + Sj_.rSj_s) = (2/m) [ns(6)ns(6') + ns(6 — 6'){ns(6) — ns(6')}] ^ Sj . (152) 



dj{sj+rdj+s + Sj^rdj-s) = "2 [{ns(6) — ns(6 — 6')}cs(6') + cs(6 — 6')cs(6)] ^ Sj . (153) 



J2 Cj{sj+rCj+s + Sj-rCj-s) = -(2/m) [{ns(6) - ns(6 - b')}ds{b') + ds{b - b')ds{b)] J2 ■ (l^^) 
i=i 3=1 
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E Cjdj(s]+r + = (Vrn) [ns'(6) + ds(&)cs(&)] cjdj 



V 



Sj(ij(sj+rCj+r + Sj-rCj-r) = (2/m) ns(6) [cs(6) + ds(6)] ^ cj 

P V 

E djCj{sj^rSj+s + Sj-rSj-s) = (2/m) ns(6) ns(6') ^ 90?^- . 



E djCj(cj+rCj+s + Cj-rCj-s) = (-2/m)ds(6)ds(6') ^ c^d^- . 

p p 
'^djCj{dj^rdj+s + dj-rdj-s) = — 2cs(6)cs(6') E'^i'^i • 
j=i j=i 

p p 

SjCj(o?j+j.Sj+s + dj-rSj-s) — (2/m)cs(6) ns(6') ^ CjO?j . 
i=i j=i 

p p 
^ Sjcij(cj+rSj+s + Cj_rSj-s) — (2/m)ds(fo) ns(fo') ^ CjO?j . 
j=i j=i 

p p 
E CjSjO?j(sj+r + Sj-r) = -(2/m)cs(6) ds(6) ^ Cjdj . 
j=i j=i 

p p 
E Cj{d^j+r + d]^r) = 2cs(6) ns(6) ^ c,d,- . 

E + cj-r) = (2/m)ds(6) ns(6) ^ c,rf,- . 



cs'(6)ns^(6) 



E^jCj-«i(4+r- + 4-r-) = -(2/m2 
+ns2(6)ds2(6) + ds^(6)cs2(6) + 3 ns2(6)cs(6)ds(6) 



E '^i'^i 
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L = 4 



p p 



+ = -(2/m)cs(6)ds(6) 



+ (2/m2)ns2(6)[ns2(6) -cs(6)ds(6)] j^s,-, 



E + = (2/^) ns^(&) E 4 

+(2/m') [cs'(6) ns2(6) + ns2(6)ds2(6) + ds'(6)cs'(6) - 3 ns'(6)cs(6)ds(6)] £ 



CjSjdj{cj+rdj+r + Cj^rdj-r) — 2 cs(6)ds(6) E *i ~ (2/"^) cs^(6) ns^(6) 
W(6)ds'(6) + ds'(6)cs'(6) - cs(6)ds(6)(cs'(6) + ds'(6) + ns'(6))j sj 



p 

2 



s|((ij+rCj+r + (ij_rCj_r) = (2/m)cs(6)ds(6) Cjdjs] 



V 

+ (2/m2) ns2(6)[ns2(6) + 3cs(6)ds(6)] ^ c^d^- 



p p 
E + dj-rCj-s) = (2/m)cs(6)ds(6') ^ Cjcij-s^^ 

j=i i=i 

r 1 P 

+(2/m2) [ns(6)ds(6)ns(6')cs(6') +cs(6)ds(6'){ns2(6) +ns2(6')}] E^i^^j 



Examples belonging to the class MI - IV 

L = : 



E dj{Sj+r + Sj-r) — 
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L = 1 : 

The generalizations of Eq. ()84|) and Eq. (jH3j) pertinent to this class are for / odd (/ < p): 

p 

X] ^j'^j+r ■ ■ ■ Cj+il-l)r = (172) 

i=i 

1 /(;-l)/2 a-l)/2 Z \ p 

-i^ n ds^(fc6) + 2(-l)('-)/^ E n ds{[n-k]b)]Y.c,. 

\ fc=l fc=l n=l,n^k ) i=l 

When / = p, the resulting identity takes the simpler form 

AKn 

j=\ ri.=l ^' j=l 



m(--^)/^nc. = n ds^(^) E^.- (173) 



E c,2(c,+. + c,_.) = (2/m)[ns(6)cs(6) - ds^(6)] ^ c, . (174) 
i=i i=i 

V V 

dj{dj-^.rCj+r + dj-rCj-r) = — 2ds(6) [cs(6) — ns(6)] ^ Cj . (175) 
p p 

Sj(Sj+rCj+r + Sj_rCj-r) = ( — 2/m)ds(6) [cs(6) — ns(6)] ^ Cj . (176) 

i=i i=i 

p p 
^ dj{cj+rdj+s + Cj.rdj-s) = -2 [{ds(6) - ds(6 - b')}cs{b') + cs{b - b')cs{b)] J2 ■ (1^7) 



J2 Sj{cj+rSj+s + Cj-rSj-s) = (2/m) [{ds(6) - ds(6 - b')} ns{b') + ns{b - b') ns(6)] ^ Cj . (178) 
i=i i=i 



p p 
^ 9(sj+^sj+, + Sj_,.Sj_,) = (2/m) [ns(6) ns(6') + ns(6 - b'){ds{b) - ds(6')}] E Cj . (179) 
j=i i=i 



Cj{dj+rdj+s + dj_rdj_,) = -2 [cs(6)cs(6') + cs(6 - b'){ds{b) - ds(6')}] E • (180) 



p p 
^ 9(cj+^,9+, + 9_,.Cj_,) = -(2/m) [ds(6)ds(6') + ds(6 - b'){ds{b) - ds(6')}] E ■ (181) 
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L = 2 : 



p 


r 1 ^ 

-(2/m) [ds^(6) + cs(&) ns(6)] Sjd^- . 


(182) 


P 

CjSjCij[Cj-(_r. + Cj_r 


= (2/m)cs(6)ns(6) ^ Sj 


dj . 


(183) 


P 


= -(2/m)ds(6) [ns(6) + cs 


p 

(6)] E Sjdj . 
j=i 


(184) 


P 


p 

= (2/m)cs(6)ds(6) ^s^- 


dj . 


(185) 


P 


p 

— (2/m)ds(6)ns(6) Sj 


dj . 


(186) 


p r 

^9d,s,(c|+, + c5_,) = -(2/m2) cs=^(6)ns'(6) 
j=i '- 

W(6)ds^(6) + ds^(6)cs=^(6) + 3ds\b) ns(6)cs(6) 


p 

J2sjdj ■ 


(187) 


^jdj{Cj-\-rCj-^s + Cj-rCj. 


_,) = -(2/m)ds(6)ds(6')E^,ci, . 


(188) 


P 

^jdjidj-^rdj+s + 


rdj-s) = -2cs(6)cs(6') 


Sjdj . 


(189) 


p 


p 

_s) = (2/m) ns(a) ns(a') ^ Sjdj . 


(190) 


P 

Cjdj{Sj-^-rCj^s + ^j-rCj- 

j=i 


p 

= -(2/m)ns(6)ds(6')IIsj^ij- • 


(191) 


P 


_,) = -(2/m)cs(6)ds(6')Es,ci, . 


(192) 
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s]d]{cj+r + Cj-r) = -2 ns(6)cs(6) J2 4 



3=1 



+ (2/m)cs(6) ns^(6) [msn^{b) + 01"^ (b) - cn(6)] J^^i ■ (1^3) 



+ (2/m2) [cs2(6) ns\b) + ns\b)ds\b) + ds'(6)cs2(6) - 3ds\b)cs{b) ns(6)] ^ . (194) 



^ CjSjdj{sj+rdj+r + Sj_rdj-r) = 2cs(6)ns(6) ^ + (2/ 
i=i i=i 



cs'(6) ns'(6) 



+ ns'(6)ds'(6) + ds'(6)cs'(6) -cs(6) ns(6)(cs'(6) + ds'(6) + ns'(6)) 



E^- (195) 



Appendix B: Examples Using Master Identities With Alter- 
nating Signs 

Note that in this case the identities are only vaUd when p is even integer and since r and p are 
coprimes hence r is necessarily odd. The letters a, a' again stand for 2rK/p and 2sK/p respectively. 

Examples belonging to the class MI - I 

L = l : 



E(-ir'^.(c.+- + c.-.) = 0, (196) 
p p 

J2{-iy~%dj+rdj+2r = -[cs^(a) + 2cs(a)cs(2a)] E(-l)^'~^^^j (197) 



E(~l)'' ^djdj+rdj+s — — [cs(a)cs(a') + cs(a)cs(a' — a) — cs(a')cs(a' — a)] E(~l)'' (198) 



Y{-'^y~^d]{dj+r + dj-r) = 2[ds(a)ns(a) +cs2(a)] ^(-1^"^^^ • (199) 
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p p 

J^i-'^y'^Cjicj+rdj+r + Cj_rdj-r) = -(2/m)cs(a) [ds(a) + ns(a)] . (200) 

j=i j=i 

p p 

Y,{-'^y~^Sj{sj+rdj+r + Sj_rdj-r) = (2/m)cs(a) [ns(a) + ds(a)] J2{-iy-^dj . (201) 
j=i j=i 



^{-iy-^dj{dj+rdj+s + dj_rdj-s) = -2[cs(a)cs(a') -cs(a-a')(cs(a) -cs(a')] ^(-l)^"^c^j- . (202) 



^(—1)"' ^dj{cj+rCj+s + Cj-rCj-s) = — (2/m) [ds(a)ds(a') — ds(a — a')(cs(a) — cs(a')] ^(—1)'' "^rfj • 
j=i j=i 

(203) 



^(—1)"' ^dj{sj+rSj+s + Sj-rSj-s) = — (2/m) [ns(a)ns(a') — ns(a — a')(cs(a) — cs(a')] ^(— l)-* ^dj . 

(204) 



P 

= (2/m)[ns(a)cs(a') - ns(a)cs(a - a') + ns(a')ns(a - a')] ^(-l)^'"^dj . (205) 



^( — 1)-' Cj(Cj+rC^j+s + Cj-rdj-s) 

i=i 

p 

= -(2/m)[ds(a)cs(a') - ds(a)cs(a - a') + ds(a')ds(a - a')] ^(-1)^-^^ . (206) 



L = 2 



X:(-l)^-M|[c,+,s,+, + Cj_rSj-r] = 2[cs'{a) - ds(a) ns(a)] f^(-l)^-ic,-5, . (207) 



^(— 1)-^ ^Sjdj{cj+rdj+r + Cj-rdj-r) = —2cs{a)[—'ns{a) + ds{a)]'^{—iy ^cjSj . (208) 

32 



Cj-rCj- 



-(2/m)ds(a)ds(a') ^(-l)^'^^CjSj- . 



p 



p 



J-sJ 



-2cs(a)cs(a')E(-l)'"'95. ■ 



p 



p 



^(-1)^ ^Cjdj(sj+rdj+s + Sj^rdj-s) = -2ns(a)cs(a') ^(-l)^ ^Cjs 

3=1 j=l 



(209) 
(210) 
(211) 
(212) 
(213) 



j=i 



L = 3 : 



^(— 1)-^ ^d'^{sj^rCj+s + Sj-rCj-s) = 2[cs(a)ds(a — a') — cs(a')ns(a — a')] ^(—1)-^ ^CjSj . (214) 

(215) 

(216) 
(217) 



^CjSjdj{dj+r + dj-r) ^'^ds{a)ns{a)J2i-'^y ^Cj^j ■ 
f:(-l)^-^.,(c|+, + c|_,) = {2/m)cs{a)ds{a)j:{-iy-'c,s, . 



3=1 
P 



3=1 

P 



E(-l)'"'c,(4+. + ^3-r) - -(2/m)cs(a)ns(a) ^{-ly-'cjs, 

3=1 3=1 



J2(-^y ^CjSjdj{sj+rCj+r + Sj-rCj-r) = (2/m)ds(a)ns(a) |^(-1)^ 

3=1 3=1 



cs^(a)ns^(a) + ds^(a)ns^a) + cs^(a)ds^(a) + ds(a)ns(a)(cs^(a) 



+ds^(a) + ns^(a)) 



3=1 



Examples belonging to the class MI - II 



L = l : 



(218) 
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While there are no ordinary identities of this class with L — 1, alternating identities abound. 

They are therefore unique and characterized by the appearance of the Jacobian zeta function. Also 
they have helped us in finding identities for the product of p sn's as well as of p en's. 



p 



Y,{-iy-%dj+r = -2cs(a) Y.i-'^y-'Zj . (219) 



p p 



T.i-^y-%^^+r = (2/m)ns(a) ^{-ly-'Z, . (220) 

j=l 3=1 



P P 



Ei-^y-'cjCj^r = -(2/m)ds(a) E(-l)'"'^. ■ (221) 



It may be noted that the above three identities are vahd for p > 4. 



f2{-'^y~^djdj+rdj+2rdj+3r = 2 [cs(a)cs(2a)cs(3a) + cs^(a)cs(2a)] ^(-l)^-^Zj . (222) 
This generalizes for any even number Z < p to: 

P (1/2 1 \ P 

Y.{-iy-%d,^r---d,^^i-,y^2{-lf'{Y^{-lf-' n cs{[n-k]a)W{-iy-'Z,. (223) 

j=l \k=l n=l,n^k ) j=l 

Similarly, for any even number I < p, the sn and cn functions satisfy the identities (p > 4) 



P / V2 i 



j=l \k=l n=l,n^k J j=l 



^■'-E(-irV.4-.---«.+{^-i).-2 X:(-l)'-' n ns{[n-k]a)\Y.{-iy-'Z,. (224) 



P i in I 



m~'-Y.{-iy-'c,c,^r ■ • ■c,+(,_i). = 2(-l)'/2 E(-l)'"' n ds([n - k]a) E(-ir'^. • 



^1/2- 

j=l \fe=l n=l,n^k J j=l 



(225) 

When / = p (p > 4) the above identities reduce to: 

n n ns^(^) E(-l)^-^^.- ■ (226) 



TT «. = I TT n«2r 

3=1 \n=l ^ / 3=1 



fp/2-l 

ds^(- 

3=1 \n=\ ^ I J=l 



m^'^ n = vi^{-ir" n ^a-^) u-^y-'z, . mi) 
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Thus these are the even product equivalents of Eqs. p46|) and fjl73|) . It is worth reminding here 
that the product did2...dp is on the other hand a constant being equal to (1 — m)^/^. 

p p 

^(-iy-^c/j[cj+^ Sj+r + Cj-r Sj-r] = -(4/m)ds(a) ns(a) J2{-iy-^Zj . (228) 

j^i-'^y-'sjicj+rdj+r + Cj.rdj-r] = -(4/m)ds(a) cs(a) Y.i-'^y^'Z, . (229) 
i=i i=i 
p p 

'^{—iy'^^Cj[sj+rdj+r + Sj-rdj-r] = — (4/m)cs(a) ns(a) '^(—ly^^Zj . (230) 
j=i j=i 



^(— !)■' ^dj[cj+r Sj+s + Cj-r Sj^s] = i4:/m)[cs{a) ds{a — a') — cs{a')ns{a — a)]'^{—iy ^Zj . (231) 
i=i i=i 



^(— !)■' Cj+s+'^j-r Cj_s] = — (4/m)[ns(a) ds(a— a')— ns(a') cs(a— a')] ^(— !)■' ^Zj. (232) 

i=i i=i 



^(—1)-' ^Cj[c?j+.r Sj+s+dj_r = — (4/m)[ds(a) ns(a— a')— ds(a') cs(a— a')] ^(— !)■' . (233) 
L = 2 : 

f](-l)^-i4(rfj+, + rfj_,) = 2ns(a)ds(a) ^^(-l)^''^^ . (234) 
i=i i=i 

E(-l)^-^ c]{c,^r + 9-.) = (2/m^)cs(a) ns(a) • (235) 

f](-l)^-i s5(s,+, + s,.r) = (2/m2)cs(a) ds(a) X:(-l)'"'ci • • (236) 



= (4/m')[ns'(a)(cs'(a) + ds'(a) + cs'(a)ds'(a)] f](-l)^-^4 . (237) 



L = 3 : 
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^ X dj\Cj-\-rSj-^r ~^ Cj—rSj—r\ 

P P 

= (12/m)cs2(a)ds(a)ns(a) -2ns(a)ds(a) J2{-iy-^SjCjdj . (238) 



^( — 1) CjSjdj[Cj+r + Cj-r] 

P P 

= (2/m) cs(a) ns(a) ^(-1)^"^ c^s^dj- - {4/m^)ds^{a) cs(a) ns(a) ^(-ly^^Z^- . (239) 



^( — 1) SjCjdj[Sj+r + Sj-r] 
3=1 

P P 

= -(2/m) cs(a) ds(a) ^(-1)^"^ c^s^rfj + {A/m^)Tis^{a) cs(a) ds(a) ^(-l)^-^Zj-. (240) 

3=1 j=i 



''i^3'^3\y3+i' ' "'J-»'J 

p p 
- 2 ds(a) ns(a) ^(-1)^'"^ CjSjdj - {A/m)cs^{a) ds(a) ns(a) ^(-1)^-^^^-. (241) 

jZ{-iy-' d,c,sAd]^, + = -2 cs2(a) i:(-l)^-^ c,.,c^,- 

+ (4/m)[ns2(a)(cs2(a) + ds'(a)) + cs^{a)ds\a)] Y,{-iy-^Zj. (242) 



i:(-i)^-^5,c,[4^,+4 



2 ds(a) ns(a 

) E(-iy"' - (12/m)cs2(a) ds(a) ns(a) ^(-l)^-^Z,-. (243) 



"3 1 
3=^ 



(2/m) cs(a) ns(a) £(-1)^"^ CjSjdj - {12 / m^)ds^ {a) cs(a) ns(a) ^^(-l)^-^^^-. (244) 
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p p 
= -(2/m) ds(a) cs(a) ^(-1)^"^ c^s^fij + {12/m^)ns^{a) ds(a) cs(a) ^(-l^-^Zj. (245) 



^(-1)^-^ rf,c,.«5[4+, + = -(S/m^) ns^(a)ds(a) cs(a) ^(-1)'"' 

r "I ^ 

3(ns^(a) + ds^(a) + cs^(a)) + cs^(a)ds^(a) 



+(4/m'^)ns^(a) ds(a) cs(a) 



^(-ly-^Z,. (246) 



j:{-iy-' d,s,cy^^, + cJJ = -{8/m') ds2(a)ns(a) cs(a) 9^.^. 



+ (4/m )ns(a) cs(a) 



3ds (a)(ns (a) + ds (a) + cs (a)) + cs (a)ns (a 



^(-l)^-iZ,. (247) 



5:(-l)^-i c,^,c?2[4^^ + 4 J = -8cs2(a)ns(a) ds(a) El-l)'"' 



j=i i=i 



+(4/m)ns(a) ds(a) 



3cs^(a)(ns^(a) + ds^(a) + cs^(a)) + ds^(a)ns^(a 



E(-1)^-^Z,. (248) 
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